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ABSTRACT 


The chameleon gravity model postulates the existence of a scalar field that couples 
with matter to mediate a fifth force. If it exists, this fifth force would influence the 
hot X-ray emitting gas Ailing the potential wells of galaxy clusters. However, it would 
not influence the clusters weak lensing signal. Therefore, by comparing X-ray and 
weak lensing profiles, one can place upper limits on the strength of a fifth force. 
This technique has been at tempted before using a single, nearby cluster (Coma, 2 : = 


0.02, Terukina et al. 2014). Here we apply the technique to the stacked profiles of 


58 clusters at higher redshifts (0.1 < 2 ; < 1.2 ), including 12 new to the literature, 
using X-ray data from the XMM Cluster Survey (XCS) and weak lensing data from 
the Canada France Hawaii Telescope Lensing Survey (CFHTLenS). Using a multi¬ 
parameter MCMC analysis, we constrain the two chameleon gravity parameters 
and 0oo)- Our fits are consistent with general relativity, not requiring a fifth force. 


In the special case of f{R) gravity (where P = a/1/6), we set an upper limit on the 
background field amplitude today of |/ro|< 6 >< 10“^ (95% CL). This is one of the 
strongest constraints to date on |/ro| on cosmological scales. We hope to improve this 
constraint in future by extending the study to hundreds of clusters using data from 
the Dark Energy Survey. 
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An accepted explanation for the accelerated expansion of the 
late-time Universe ( Riess et al.|1998' Perlmutter et al.|1999 ) 
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is to modify the Einstein equation, either by adding a com¬ 
ponent to the energy-momentum tensor via dark energy, or 


to the Einstein tensor via a modification to gravity (Milgrom 
Clifton et al.|2012 ). The latter often involves the intro- 


1983 


duction of a scalar held coupled to the matter components 
of the Universe, giving rise to a hfth force of the same order 
of magnitude as gravity ( Jain et ar]|2013 ). Through a vari¬ 
ety of experiments and astronomical observations, this hfth 
force has been demonstrated to be negligible at terrestrial 
and solar system densities ( Wagner et al.||2012 ). Therefore, 
if a hfth force does exist it must be it must be suppressed, 
or “screened”, in high density regions and only take ehect 
in low density regions. 

One model with such a screening is the chameleon mech¬ 


anism (Khoury & Weltman 2004). In this approach, the 
scalar held coupling strength is sensitive to the depth of the 
local gravitational potential. In regions with a large poten¬ 
tial well this screening suppresses the hfth force and gravity 
behaves as predicted by general relativity. However when the 
potential becomes small, the hfth force is unsuppressed and 
gravity becomes “modihed” compared to general relativity 
( Lombriser||20I4 ). 

By dehnition, the chameleon held satishes 




( 1 ) 


( Khoury Weltman||2004 ), where V is the potential of the 
scalar held, /3 is the coupling between matter and the scalar 
held, 0 gives the position dependent screening efficiency, Mpi 
is the Planck mass and p is the matter density. This leads 
to the chameleon hfth force of 




( 2 ) 


There is a particular set of gravity models, known as 
f{R) models ( [Buchdahl 1970) which exhibit a chameleon, 
where the strength of the hfth force (parameterised by jS in 
Equation has a hxed value (3 — a/I/ 6. This force arises 
from adding a scalar function f{R) to the Ricci scalar in the 
Einstein-Hilbert action ( Capozziello|2QQ^|Nojiri Odintsov 
2003). These models can reproduce observed late time ac¬ 


celeration of the Universe whilst still suppressing the hfth 
force in high-density environments, such as the Solar Sys¬ 
tem ( Chiba et al.|2007 ). These f{R) models possess an extra 
scalar degree of freedom, /r = d//dR, where the value at 
the current epoch is |/ro| ( [Sotiriou &: EaraonT||20I0 ). Then 
f{R) gravity can be related to 0oo, (0 in Equation at in- 
hnity) via the relation ( [Joyce et al.||20I4| ) 




(3) 


Hu & Sawicki (2007) provide theoretical arguments 


showing that for general relativity to be preserved at parsec 
scales within the Solar System , then |/ro|< I0“®. At kilo- 
parsec scales Jain et ^[2013 constrained |/ro|< 5 x 10“^ 
in dwarf galaxies. On megaparsec and larger scales, [Raverij 
et al. ( |20I4| ) used the Cosmic Microwave Background to 


/ro|< 10 Also on large scales Rapetti et al. 
(2011), [Ferraro et ~ (2011), Cataneo et al.| ( |2QT4 ) used the 


abundance of galaxy clusters to constrain |/ro|, e.g. Cata-j 
jneo et al. (2014) measured (under the assumption of n = I), 

|/ro|<2.6x 10-^ 


In this paper, we also use clusters of galaxies to con¬ 
strain |/ro| on megaparsec scales. However, unlike 
et al.| < |2011| ), [Eerraro et al. ( |20II| , [Cataneo et al. 


Rapetti 


(2014), 


we use cluster prohles, rather than abundances to do so. 
The hypothesis is that a fifth force would be screened in the 
dense cluster cores, but not in the rarefied cluster outskirts 


( Burikham Panpanich[[20I2 Lombriser et al.[|20I2 ). The 
majority of baryonic matter in a cluster is ionised gas that 
has been pressure-heated to temperatures in excess of lO^K 
( Gursky et al.|I9'7T Loewenstein|2004 ), leading to the emis¬ 
sion of X-rays via thermal bremsstrahlung radiation (Jones 
&: Eorman|I978 Sarazin|2009 ). The gas can also be observed 
indirectly through its influence on the cosmic background 
radiation, via the so-called Sunyaev-Zel’dovich (SZ) effect 
( Sunyaev Zeldovich[[I98Q ). 

By measuring the properties of this X-ray gas we are 
able to infer, under the assumption of hydrostatic equilib¬ 
rium, the cluster mass and density from its X-ray surface 
brightness or SZ effect profiles ( Reiprich Bohr inger| 2002 


Kettula et al.|2QT4 ). In a chameleon gravity model, the intr¬ 
acluster gas would feel the fifth force in addition to gravity 
in the cluster outskirts, i.e. the gas will be slightly more 
compact and the temperature boosted ( Arnold et al.|2014 ), 
compared to the influence of general relativity alone. 

By contrast, weak gravitational lensing is dependent 
only upon the gravitational deflection of light by matter 
along the line of sight, therefore providing a technique to 
measure the underlying mass distribution without assum¬ 
ing hydrostatic equilibrium. Crucially for this study, the 
fifth force would not modify the deflection of light through 
the cluster (compared to general relativity) because the 
scalar chameleon field is coupled to the trace of the energy- 
momentum tensor ( [Hui et al.[2009[ ). Therefore, we can search 
for evidence of a fifth force by comparing the X-ray sur¬ 
face brightness, and/or SZ effect, profiles of clusters with 
their gravitational lensing shear profiles (Ostriker & Vish- 
I986[ [Terukina Yamamoto[[20I2 ). 


Terukina et al. 


^ 20I4[ ) used this approach to constrain 


f{R) gravity models using a combination of lensing shear. 
X-ray surface brightness. X-ray temperature and Sunyaev- 
Zel’dovich profiles for the Coma cluster (a massive cluster at 
z = 0.02). Combining these measurements, they performed 
an MCMC analysis of the parameter space describing the 
cluster profiles in the modified gravity regime. Under the 
assumption of hydrostatic equilibrium, they obtained con¬ 
straints of |/ro|< 6 X 10“^. They also examined the assump¬ 
tion of hydrostatic equilibrium, and concluded that any con¬ 
tribution of non-thermal pressure was small compared to the 
reconstructed mass. 

The Coma cluster is at low redshift meaning its weak 
lensing shear signal is low. Moreover it is known to have non- 
spherical geometry (Eitchett & Webster I987[ Briel et al. 
1992 Colless I)unn[[I996 ). These factors motivate us to 
apply the Terukina et al. (2014) method to many more clus¬ 


ters at higher redshifts, allowing for a higher signal to noise 
weak lensing shear profile and an averaging out of non- 
spherical cluster shapes. We do this by comparing stacked 
X-ray surface brightness and shear profiles of 58 X-ray se¬ 
lected clusters. We utilise high quality weak lensing data 
from the Canada Erance Hawaii Telescope Lensing Survey 
(CEHTLenS, Hey mans et al.[[20I2 Erben et al.[[20I3 ), and 
X-ray observations from the XMM Cluster Survey (XCS, 
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Romer et al.|200T]|Lloyd-Davies et al.|20lT]|Mehrtens et aL] 


2012| . We also investigate the Temkina et al.| ( |2014| con- 
clusion that deviations from hydrostatic equilibrium do not 
invalidate the chameleon gravity test. 

In Section!^ we review the underlying theoretical back¬ 
ground. In Section we describe the development of the 
cluster sample used in the analysis, and the Markov chain 
Monte Carlo (MCMC) methods used to simultaneously ht 
the X-ray surface brightness and weak leasing prohles. In 
Section we discuss our results and the implications of our 
results in the framework of f{R) gravity models. In Sec¬ 
tion [5] we discuss the influence of cluster environment and 
of our assumption of hydrostatic equilibrium. In Section 
we present our conclusions. Throughout this paper we use 
a 95% conhdence level when quoting upper limits, adopt a 
cosmology with Qm = 0.27, Qa = 0.73 and Hq — 70 km s“^ 
Mpc“^. 


2 THEORETICAL BACKGROUND 

In this study, we adopt the Navarro-Frenk-White (NFW; 
Navarro et al. 1996) model for the dark matter halo mass 


distribution: 


p{r) = 


Pc^c 


— (1 + —) 


(4) 


where r here and throughout is the radial distance from 
the halo centre, pc = 3H‘^{z)/87tG is the critical density at 
a given redshift, H{z) is the Hubble parameter at a given 
redshift, G is Newton’s Gravitational Constant, 5c is the 
characteristic overdensity, given by 

200 


(5c = 


(5) 


3 ln(l + c) — c/(l + c) ’ 

where c is a dimensionless concentration parameter and r. 
is the scale radius given by 


^ 1 

C\A-Kpc5c) 


( 6 ) 


where M 200 is the mass enclosed by r 2 oo, the radius at which 
the dark matter halos average density is two hundred times 
the critical density. 


M(< r 2 oo) = 47r5c/5c?'a ^ln(l + c) - 


(7) 


The NFW prohle described in Equation is well sup¬ 
ported by N-body simulations of ACDM, but it is not im¬ 
mediately obvious that this prohle would pertain to clus¬ 
ter prohles in the f{R) regime. However it has been shown 
(Lombriser et al. ] 2012| [Moran et al.|[MT^ that the NFW 
prohle is able to provide hts to both modihed gravity and 
concordance cosmology that are equally good, shari ng the 
same It should be noted that the simulations in Lom¬ 


briser et al. 


(2012) were generated using a hxed (3 — y/l/6. 


as opposed to the general chameleon gravity model inves¬ 
tigated here. However as we are probing a /3 range around 
this value, we expect any modihcations to the prohles to be 
similar, suggesting the suitability of the NFW prohle. Fur¬ 
ther checks using hydrodynamical simulations of modihed 
gravity models would allow this assumption to be verihed. 


We adopt the Temkina et al. (2014) approach describ¬ 


ing the chameleon mechanism using three parameters. The 


hrst of these, /3, is the coupling between matter and the 
scalar held (see Equation [^. The second, 0oo, describes the 
position dependent screening efficiency. The third, rcrit, is a 
critical radius, i.e. the distance from the dark matter halo 


centre at which the screening mechanism takes effect (Teruk- 
ina Yamamot^|2012 ), 


^crit — 


PPstI 

Mpi0o 


(8) 


where ps is the density at this radius. 

[Temkina Yamamoto (2012) showed the hydrostatic 
equilibrium equation in the presence of a hfth force (Equa¬ 
tion]^ is: 

1 dPgas(r) _ GM{< r) ^ 


Pgas(r’) dr 




where pgas is the gas density, M the total mass within a 
radius r, and Pgas is the electron pressure. 

In an ideal cluster, i.e. one that is isolated, isother¬ 
mal, and spherical, this total pressure is felt by the elec¬ 
trons and ions in the ionised intracluster plasma, so that 
Pgas = riekT, where rie is the electron number density, and 



nr 

Pe{r) = Pe,o + pmp / ne(r) 

Jo 

( GM{<r) _ jJ d<p{r) \ 

V r2 Mpi dr ) ' 


( 10 ) 


where Pe,o is the electron gas pressure at r = 0, given by 
Pe,o = ne,o^T and ne,o = 5no/(2 + p) and M(< r), the 
halo mass. The integral of Equation can be re-expressed 
in terms of a projected X-ray surface brightness SbC) us¬ 
ing the temperature and electron density dependent cooling 


function (see Section 3.2), 


SB{r±) = 


1 


47r(l + z)^ 


j (% 


rl+z‘^]\o{T^^)dz, ( 11 ) 


where r± is the projected distance from the cluster cen¬ 
tre and z the cluster redshift. This is the expression we 
ht to when comparing stacked X-ray cluster prohles to the 
chameleon model (Section [3.5[ ). 

The expression used to ht the weak lensing shear prohles 
(under the assumption of an underlying NEW prohle) for 
comparison is given in [Wright &: Brainerd|2QQ0l 

To recap, our method makes the following assumptions: 
that modihcations to general relativity include a chameleon 
screening mechanism and can be described by Equation H 
that dark matter halos follow an NEW prohle (Equation 
that a hfth force can be included in the hydrostatic equi¬ 
librium expression according to Equation [^ that clusters 
of galaxies are isolated, isothermal, and spherical (which in 
turn implies that the clusters are in hydrostatic equilibrium, 
have an electron number density that follows a beta-model, 
and their X-ray emission can be predicted from a thermal 
cooling function); and that the weak lensing shear prohles 
of clusters are given in [Wright Brainerd|200^ We discuss 
the impact of some of these assumptions in Section]^ 


1 


The beta in this model is not the same as the /3 in 


Equation]^ 
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3 METHODS 

3.1 Compiling the X-ray Cluster Sample 

In this paper we used public weak leasing data (galaxy ellip- 
ticities and photometric redshifts) provided by the Canada 
France Hawaii Leasing Survey (CFHTLenS, |Heymans et al. 


2012). The CFHTLenS covers 154 sq. deg with high quality 


shape measurements. The galaxy ellipticities were generated 
by the CFHTLenS team using the THELI ( Erben et al.|2013 ) 
and lensfit ( [Miller et al.|2013 ) routines. Photometric red- 
shifts were produced using PSF-matched photometry to an 
accuracy of 0.04(1 + z) with a 4% catastrophic outlier rate 
( Hildebrandt et aL]|2012 ). 

We also used public X-ray data taken from the XMM- 
Newton archive and have collated a sample of X-ray clusters 
in the CFHTLenS region using pipelines developed for the 
XMM Cluster Survey (XCS, [l7loyd- Davies et al.|2011 ). First 
we determined which of the XMM observations overlapped 
with the CFHTLenS fields. We then used the XCS pipelines 
to carry out the following tasks in an automated manner: 
cleaning the event lists of background flares; creating de¬ 
tector and exposure images; producing duplicate free source 
lists; and identifying extended X-ray sources. A total of 348 
extended XMM sources, with more than 100 background 
subtracted photon counts, were located in the CFHTLenS 
helds, although 44 were close to the edge of the XMM held 


of view and were not considered further (please see Lloyd- 


[Davies et al.||20lT] for the relevant, XCS specihc, dehnition 
of source counts). 

The majority of these sources were not included in the 


XCS hrst data release (XCS-DRl, Mehrtens et al. 2012). 
This meant that candidate identihcation needed to be car¬ 
ried out before the sources could be used in our study. This 
process is non trivial: as shown in [Mehrtens et al.|[2012[ a 
large fraction of XCS extended sources (especially those with 
fewer than 300 counts) are either hard to conhrm as clusters 
- because the available imaging is not deep enough - or are 
associated with other types of X-ray source. Therefore, for 
this paper, we have taken a conservative approach and only 
included XMM extended sources in our study if they corre¬ 
spond to an over density of galaxies in false colour images 
produced using the CFHTLenS cutout servic^ One hun¬ 
dred and eighty six sources were excluded from the study as 
a result. These were excluded for several different reasons: 
there was no optical data as the cluster sat in a masked re¬ 
gion of the CFHTLenS footprint; there was a bright star or 
galaxy lying close to the cluster centre that was obscuring 
it; or the optical image resembled an AGN rather than a 
cluster. The coordinates of the remaining one hundred and 
nineteen can be found in Table IbTI 

As our analysis required information about the distance 
to the cluster, a further 37 sources were excluded from the 
study because redshifts were not available at the time of 
writing. These are flagged with a 2 in Table [BT] The major¬ 
ity (63 of 82) of the redshifts we used came from the new 
Gaussian mixture model redshift estimator described in de¬ 


tail in IHood X Mannll2015l We also use d 18 redshifts taken 
from NEe[^ and 3 from [Ford et ^ ( [2014 ). 

We judged these remaining 82 XMM extended sources 
in the CFHTLenS region to be conhrmed clusters and ran 
them through the XSPEC based XCS spectral pipeline. We 
determined X-ray temperatures when the signal to noise was 
sufficient. This produced X-ray temperatures of 58 of these 
clusters which form our final sample, including 12 clusters 
new to the literature, the other 23 clusters were excluded 
from the analysis and are flagged with a 3 in Table M 
The details of this pipeline can be found in [ Lloyd- Davies[ 
[et al.[[20TT] These 58 clusters with measured temperatures 
span the redshift range 0.1 < < 1.2 (median 2 : = 0.33) and 

temperature range 0.2 < Tx < 8 keV (median Tx = 2.3 keV). 
A selection of these new to the literature clusters, along with 
several clusters that were optically confirmed but excluded 
due to a lack of redshift, are shown in Eigure[Cl 


3.2 Making stacked X-ray Surface Brightness 
Profiles 

Our analysis involves stacking multiple different XMM ob¬ 
servations of our 58 clusters, in order to build up signal 
to noise in the outer parts of the ensemble cluster prohle. 
This process needs to account for the following complexities: 
Most of the 58 clusters were covered by more than one XMM 
observation. Each of these observations has different back¬ 
ground properties and flare corrected exposure times. The 
X-ray telescope comprises of three cameras that operate si¬ 
multaneously (mosl, mos2, pn), so most XMM observations 
comprise of three separate images with different, energy de¬ 
pendant sensitivities. The clusters all have different energy 
spectra, because, even if one ignores non thermal processes, 
they have different X-ray temperatures, redshifts, and line of 
sight absorbing column densities. Therefore, for each cluster, 
we have to calculate, using XSPEC, camera specihc count 
rate to luminosity conversion factors for each XMM obser¬ 
vation that it falls in. We then, for a given cluster, take 
the photon count images generated by the XCS pipeline, di¬ 
vide these by the respective exposure map, and multiply by 
the cluster dependent conversion factor. This allows us to 
combine all the images for that cluster in a self consistent 
manner. 

To produce a single stack, we hrst re-scaled the 58 com¬ 
bined images of individual clusters to a standard projected 
size. Eor this we estimated M 500 , the mass enclosed within 
a sphere at which the average density is 500 times the criti¬ 
cal density, using the prescription described in [Sahlen et al.[ 
(2009). A conversion between M 500 and M 200 was made 


following the formulae derived in Hu X Kravtsov (2003) 


where we assume c = 5. This is an accurate description of 

and is 


the typical density prohles in clusters (Arnaud 


consistent with the hndings of Kettula et al. ( 2014[) in the 


2005 


CEHTLenS region. Using the M 200 values we calculated the 
radius at which the average density is two hundred times the 
critical density, r 2 oo, following the method in [Croston et al.[ 


^ http://www.cadc-ccda.hia-iha.nrc- 
cnrc.gc.ca/community/CFHTLens/cutout.html 


^ The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech¬ 
nology, under contract with the National Aeronautics and Space 
Administration. 
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(2008). The 58 stacked images could then be rescaled using 


linear interpolation to a common 500 by 500 pixel format, 
so that they each had an r 2 oo radius of 125 pixels. Each of 
these 500 by 500 images were centred on the source centroid 
as determined by XCS. 

We re-scaled the individual cluster images by the over¬ 
all amplitude of their X-ray surface brightness, as adding 
clusters over a range of different masses and luminosities 
would result in significant off-diagonal elements in the co- 
variance matrix of the final stacked prohle. Therefore, we 
calculate the mean value of the X-ray surface brightness pro¬ 
hle for each cluster, and re-scale individual cluster surface 
brightness maps by this value (we found that using the me¬ 
dian value instead of the mean gave similar results). A final 
stacked surface brightness map of the 58 individual clusters 
is then produced by taking the mean value for each pixel 
across all these maps. This re-scaling of the amplitudes is 
permitted as our constraints on modihed gravity parameters 
focus on the shape of the cluster prohles; we marginalise over 
the amplitudes of the stacked X-ray surface brightness pro¬ 
hles in Section 4. The error covariance matrix of the stacked 
prohle was then measured directly. 


3.3 Making Stacked Weak Leasing Profiles 

We outline here the procedure used to obtain the stacked 
cluster shear prohle, 7 t, using source galaxies from 
CFHTLenS. The CFHTLenS catalogue provides measure¬ 
ments of both ellipticity components (ei and 62 ), as well as 
photometric redshifts for each source galaxy. Before shears 
can be derived from these quantities, small multiplicative 
and additive corrections {m and C 2 ) must be applied, derived 
from the dataset. We calculate C 2 and m for each galaxy as a 
function of size and signal to noise (using Equation 17 and 19 
. Hey mans et al.|2012 ). Each galaxy was weighted with the 


CFHTLenS catalogue WEIGHT parameter and calibrated 
by 


6int,i 


Oi ^2,i 
1 + m 


( 12 ) 


where C 2 was applied on a galaxy by galaxy basis and m is a 
summation of 1 + m for each galaxy, applied as an ensemble 
average to each radial bin (discussed below). 

We have an effective galaxy density, rieff, (Heymans 


et al. 2012) of 12 galaxies per square arcminute. In order 


to minimise the contamination between the lensed galaxies 
and the cluster members, we only use source galaxies with a 
photometric redshift greater than ^^duster + 0 . 2 . Our redshift 
cut is made so that there is negligible contamination be¬ 
tween cluster and source galaxies. The photo-z cut does not 
require a redshift dependence as the photo-z errors of the 
source galaxies in CFHTLenS are approximately flat close 
to the redshift of our clusters ( Hildebrandt et al.|2012 ). 

For each galaxy we calculate the tangential and cross 
shears ( 7 t, 7 x) as a function of their position relative to the 
cluster position, via the angle 0 between the cluster and 
galaxy from a baseline of zero declination. The tangential 
shear measured around each XCS determined cluster cen¬ 
troid was binned into 24 equal spaced logarithmic annuli 


that we previously scaled the X-ray profiles in Section |3.2| 
for consistency. 

Finally, in order to improve the signal to noise of the 
tangential profiles, the 58 individual cluster profiles were 
stacked. This was achieved by summing the profiles of each 
cluster and calculating an average shear in each bin across 


all clusters (McKay et al. 2001 Sheldon et al. 2009). The 


error covariance matrix was then directly measured for our 
stacked profile. Due to the large uncertainty in the central 
bin, driven by the low number density of galaxies, we exclude 
the central 0.1 x r 2 oo- 

We perform consistency and null tests upon the 
CFHTLenS shape data to ensure our recovered profiles are 
unbiased and not artefacts of the data. Figure shows the 
tangential signal (solid blue) and the cross shear (dashed 
red) around the stacked clusters. The tangential shear sig¬ 
nal has a detection significance of > 30cr while the cross 
shear signal is consistent with zero at all radii. 

Figure shows the tangential shear (solid blue) and 
cross shear (dashed red) around 58 random stacked positions 
within the overlap of the CFHTLenS region and the XCS 
footprint. The measurements in both these cases were found 
to be consistent with zero on all scales. 

For Figure we show the tangential shear around the 
stacked clusters after we have split the source galaxies into 
three bins based upon their signal to noise ratio, S/N < 20, 
20 < S/N < 40, and S/N > 40, with similar redshift distri¬ 
butions (median redshifts of 0.85, 0.82, 0.79 respectively). 
We find that the three measurements are consistent with 
each other as expected. 

Finally Figure P shows the tangential shear around 
the stacked clusters with the source galaxies cut into three 
bins based upon their photometric redshift, z < 0 . 6 , 0.6 < 
z < 0.8 and z > 0.8. At higher redshifts there are a smaller 
fraction of cluster galaxies and galaxies in front of the clus¬ 
ters, and the weak leasing signal grows with redshift. We 
see these effects as our measured signal is strongest in the 
high redshift bin. We therefore conclude we are detecting a 
genuine weak leasing signal. 


out to a distance of 10 x r 2 oo (calculated in Section 3.2). We 


then scaled the values in each of these bins in the same way 


3.4 Binning in X-ray Temperature 

To generate tighter constraints upon the modified gravity 
parameters we split our data set into two separate mass bins 
to reduce errors caused by mixing clusters of varying sizes 
and masses. We find doing so improves our constraints on 
the modihed gravity parameters compared to using a single 
bin. We cut at an X-ray temperature of T = 2.5keV, to give 
two bins of mass with equal errors on their stacked prohles. 
We note that this temperature cut approximately cuts our 
sample into galaxy clusters and galaxy groups ( [Stott et al. 
2012). Our low temperature bin (T < 2.5keV) has a median 
redshift of z = 0.32 and is hagged with a 0 in Table |B1| 
while the other (with T > 2.5keV) has a median redshift of 
z = 0.34 and a hag of 1. We repeated the analyses with three 
and four temperature bins and found no improvement in the 
constraints on the modihed gravity parameters. Therefore 
to aid with computation, we complete our analysis with the 
simplest two bin case. 
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Figure 1. Tests around the 58 CFHTLenS stacked cluster; details are provided in the text, la: Tangential and cross shear, lb: Tangential 
and cross shear around 58 stacked random points. Ic: Tangential shear for three different signal to noise bins. Id: Tangential shear for 
three different redshift cuts. 


3.5 MCMC Analysis 


4 RESULTS 


We use MCMC (Gilks et al. 1996) to fit models to our 


stacked profiles. We allow all parameters that depend upon 
the cluster properties to vary for each temperature bin. This 
leads to a total of fourteen free parameters for the four 
stacked profiles (our measured weak lensing and X-ray pro¬ 
files in two temperature bins) used to constrain modified 
gravity. Four of these were used to model the weak lensing 
mass (defined in Equations I35§- We introduce the nota¬ 
tion I, II to indicate the temperature bins T < 2.5, T > 2.5 
respectively so c^, M^qq and M^qq are the concentration 
and mass parameters for each temperature bin respectively. 

We modelled the X-ray surface brightness, using the 
method prescribed in Section by defining, for both tem¬ 
perature bins, the electron number density (itself dependent 
upon no, no^ rl and rY), and the normalisation of 

the gas temperature Tq and Tq^. We reconfigure the param¬ 
eters as = /3/(l + /3) and 0 oo ,2 = 1 — exp(—0oo/lO“^Mpi) 
to span the parameter range of /3 and 0oo in the interval 
[0,1]. To obtain the cooling function (used in Equation EH’ 
we used the XSPEC software ( Arnaud|1996| and utilise the 
APEC model ( Smith et al.||20Ql ) over a range of 0.5keV to 
2keV, i.e. the same energy range as our observations from 
XMM. This model has as inputs the gas temperature, the 
cluster redshift, the cluster metallicity and a normalisation, 
and provides the X-ray cluster flux. We adopt a metallicity 
Z = 0.3^0 ( Sato et 'M^|2011 ) throughout. Using this model 
we generate fluxes for a range of temperatures which are 
interpolated for use in our chameleon gravity model. 

The chameleon parameters ^2 and 0oo,2 are the same 
across the two bins, as the modifications to gravity should 
be independent of the cluster’s mass. 

We performed an MCMC analysis using the emcee 


code (Eoreman-Mackey et al. 2013), which implements a 


Metropolis-Hastings algorithm (MacKay 2003). We mini¬ 


mized the goodness of fit using a x statistic derived from 
joint fitting of both models (see Appendix 0 - 

Our MCMC run was a parallelised implementation us¬ 
ing 128 walkers with 10000 time steps. We removed the first 
2000 iterations as a “burn in” phase. 


In Eigurej^we show our measured X'ray and weak lensing 
profiles for both X-ray temperature bins. Our X-ray surface 
brightness profiles have been measured out to 1.2 x r 2 oo with 
high signal to noise. Likewise for our two weak lensing pro¬ 
files we have recovered a shear signal out to 10 x r 2 oo with 
high signal to noise. Also shown in Eigurej^are our best fit 
models for t he each profile using the parameters outline d in 


Section |3.5| and minimising as described in Equation Al 


We show the 2D contours for constraints on model parame¬ 
ters in Eigure[Dl 

In Eigurej^we show the 2D constraints for ^2 and 0oo,2. 
To generate our constraints we have marginalised over the 
measured likelihoods of the nuisance parameters (those that 
are not ^2 and 0oo,2). We are able to do so as we are in¬ 
sensitive to the overall amplitude of our profiles, only the 
profiles shape matters for our constraints. In Eigure we 
also show the dashed (dash-dot) line the 95% (99%) con¬ 
fidence limit excluded region from [Terukina et al.| ( |2014| ) . 
The constraints are tighter from this work on larger values 


of 13 than in Terukina et al. (2014), whilst the constraints 


on smaller values of /3 are looser. As the profiles presented 
in this work extend further from the cluster than the Coma 
profile, we probe further outside the critical radius, Vc and 
are able to better constrain large values of (3. However, as 
the errors on the X-ray profiles (and the lack of available 
SZ data) used in this work are larger than those measured 
in [Terukina et alT (2014), we are less able to differentiate 
a chameleon profile from a GR one at lower values of f3, 
leading to less constraining power. 

The shape of the contours in Eigurej^can be understood 
by considering the meaning of the parameters used in defin¬ 
ing chameleon gravity. Recall that /3 dictates the strength 
of the fifth force and <^00 is the effectiveness of the screen¬ 
ing mechanism. Therefore at low values of /3, the fifth force 
causes a deviation to the profile which is too small to be dis¬ 
tinguished from GR given the observational errors. Likewise 
as GR gravity is recovered outside the critical radius rcrit, 
this sets an upper limit on /3/0oo- As /3 increases, a lower 
value for <^00 is required to keep rcrit within the cluster, giv¬ 
ing rise to the triangular shape of the excluded region. 
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Figure 2. X-ray surface brightness profiles (left) and weak lensing (right) for the two bins of X-ray temperature: T < 2.5keV (top) 
and T > 2.5keV (bottom), against radial distance normalised by r 2 oo, the radius at which the density is two hundred times the critical 
density. We choose to show the modified gravity profiles with the highest likelihood parameters, Tq = 12.6 keV, = 2.0 X 10“^cm“^, 
b\ = -0.42, = 0.06 Mpc, = 12.2 x lO^^M©, = 3.5, Tq” = 7.8 keV, = 4.9 x lO-^cm-^, = -0.89, = 0.05 Mpc, 

MIIq = 13.7 X lO^^M©, c” = 3.8, (3 = 2,000 = 2.1 x lO-^Mpp 


4.1 Implications for f{R) Gravity 

Our constraints have implications for f{R) gravity models, 
which contain a chameleon mechanism for which j3 — y/l/6 
( Starobinsky|2007| (shown as the vertical line in Figure]^. 

From Figure we estimate an upper bound on f{R) 
gravity of 0oo < 5.8 x 10“^Mpi at 95% confidence limit, 
and therefore using Equation]^ /r(^ = 0.33) < 4.7 x 10“^ 


at 95% confidence limit (where z = 0.33 is our cluster sam¬ 
ples median redshift). The time-evolution of the background 
/r(z) for a Hu-Sawicki follows dLi et al.||20l^ , 

/r( 2) = |/Ro|-[(l+3nA)/(nM(l + a)"+4nA)]"+\ (13) 
n 

where n is a free parameter of the model. At high redshifts, 
the background energy density is higher, therefore /r(^) is 
smaller and the screening is more efficient. So /r(^) de- 
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Figure 3. The 95% (light grey region) and the 99% confidence limit (mid grey region) constraints for the chameleon model parameters 
renormalised between [0,1], (32 = (3/{l3- (3) and (f)oo,2 = 1 — exp(—^oo/lO^'^Mpi) obtained from the MCMC analysis of our combination 
of weak leasing and X-ray surface brightness for our two cluster stacks. Above the dashed (dash-dot) line is the 95% (99%) confidence 
limit excluded region from 


Terukina et al. 


(2014). The vertical line is at ^ = ^1/6, showing our constraints for f{R) gravity models. 


creases by 22% from our median redshift {z = 0.33) to z = 0, 
when n = 1, and our constraint at z = 0 is |/ro|< 6 x 10“^ 
at 95% confidence limit. Considering a Hu-Sawicki model 
with n = 3, our constraint becomes |/ro|< 2 x 10“^ at 95% 
confidence limit. Our results are comparable to the results 
for the Coma cluster reported in |Terukina et al. (2014) of 

|/ro|<6x 10-®. 


5 DISCUSSION 

In this section we discuss the influence of local overdensities 
upon our cluster sample. We also question the validity of 
the assumptions we have made while constraining chameleon 
gravity, primarily the assumption that our cluster stack is 
in hydrostatic equilibrium. 


5.1 Influence of Cluster Environment 


In addition to self screening, a cluster may be screened by 
nearby clusters and therefore still show no evidence of mod¬ 
ified gravity, even in its outskirts. To check whether this 
was expected for any of our clusters we estimated the D 


parameter detailed in Zhao et al. (2011), a parametrisation 


of the separation between a given cluster and the nearest 
larger cluster, scaled by the given cluster’s r 2 oo- We describe 
clusters with log^Q > 1 as “isolated” and clusters with 
logj^Q D < 1 as living in dense environments, and therefore 
screened. As our X-ray clusters are an incomplete set of all 
clusters in our area, we looked at overdensities in the galaxy 



Figure 4. The minimum D parameter for each cluster against 
X-ray temperature, where log^^Q H is a measure of the distance 
between a cluster and the nearest overdensity in the top 30% 
(10%) of overdensity values, shown as a red circle (blue cross). 
The shaded region contains clusters with potential screening from 
neighbouring overdensities. The majority of the clusters are in an 
isolated region. 


© 0000 RAS, MNRAS 000 , 000-000 






















Chameleon gravity with clusters 9 


density field as a proxy for nearby clusters. We binned the 
galaxies in the CFHTLenS catalogue into 3-D pixels of vol¬ 
ume IMpc^ in area, and 0.01 in redshift. Figure]^ shows X- 
ray temperature against log^^Q D, where we have calculated 
log^o ^ values between each cluster and overdensity and se¬ 
lected the smallest log^^Q D as a measure of environment. It 
is seen that only 7% (2%) of our clusters are found to be 
near (log^Q D < 1) the most overdense 30% (10%) of the 3- 
D pixels. We therefore conclude that our sample appears to 
be largely environmentally unscreened by nearby clusters, 
and therefore will apply our analysis to the full cluster sam¬ 
ple. We note that it is possible that clusters outside the edge 
of the CFHTLenS observations could screen at most 6% of 
our sample, which lie within log^Q D = 1 of the edge. 


5.2 Assumption of Hydrostatic Equilibrium 

Even in the absence of a fifth force, the interpretation of ap¬ 
parent differences in cluster mass profiles derived from X-ray 
or Sunyaev-Zel’dovich (SZ) observations and lensing mea¬ 
surements is complicated by both astrophysical processes in 
clusters, such as gas clumping in the cluster outskirts, and 
systematic errors in the measurements themselves. This has 
led to uncertainty in mass calibration being the dominant 
source of error on cosmological constraints derived from SZ 


cluster catalogues (e.g., Hasselfield et al. 2013 Reichardt 
et al.|20l3 Planck Collaboration et al.|2014 ). The absolute 


cluster mass scale is affected by uncertainty in the effects of 
feedback from active galactic nuclei, and non-thermal pro¬ 


cesses such as bulk motions, on the cluster gas (e.g., Nagai 
et al.j2007 ). Instrumental calibration uncertainties may also 


play a role (e.g., Schellenberger et al. (2015), Israel et al. 


2015). Lensing measurements, which are affected by differ¬ 


ent systematics, are being used to quantify any bias in the 
absolute mass scale, but at present, samples are small, and 
there is some disagreement (e.g., von der Linden et al.|20f4 
Hoekstra et al.||2015 ). 


In this work, we have investigated one of these issues: 
the impact of non-thermal pressure on our conclusions about 
chameleon gravity (whilst maintaining the simplifying as¬ 
sumptions of spherical symmetry). We plan to investigate 
the other issues, using hydrodynamic simulations, in future 
publications. The thermal mass of a cluster is defined by 
the gas pressure, density and temperature, which we infer 
from the X-ray surface brightness. We follow the paramet¬ 


ric fits described in Terukina et al. (2014) to reconstruct 


the stacked cluster temperature profile and electron number 
densities from the profile parameters fit for by our MCMC. 
We infer from X-ray observations. 


_ -/cTgasr f d\nTie , dlnTgj 
-‘^thermal — 


fimpG \ d\nr d\nr 


(14) 


where k is the Boltzmann constant, mp is the proton mass. 
According to the hydrodynamical simulations in | Shaw et al.| 
(2010), the non-thermal pressure can be modelled as a 


function of the total pressure, such that Pnon-thermai(^) = 
^(r)Ptotai(r), where 


S(r) = a„t(l + z)'5”‘ (^)""‘ ( 
V^’500/ V 


M 200 


3 X 1014 Mq 


(15) 


with Qiit, /dnt, Tint and um are constants determined from 16 
simulated clusters, with a mass range between 0.35 — 9.02 x 
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lO^^M© at 2 : = 0 ( Lau et ^|2009| . We adopt their best fit 
values of /3nt, n-nt, n-M = 0.5, 0.8, 0.2 respectively. In order to 
test the robustness of our assumptions we select a = 0.3, 
which was the most extreme value found in the 16 clusters 
in their analysis. The extra mass component which would 
be inferred from X-rays due to such non-thermal pressure 
would be 




d 


non-thermal — 


Gpgns dr \1- g{r) 


9(t 


skTa, 


(16) 


where r is the radial distance, g{r) is defined in Equation 15 
and pgas,rigas and Tgas are the gas density, number density 
and temperature respectively. 

In Eigure we show our mass profiles for 0.3 Mpc 
< r± < 2 Mpc for the lensing mass and X-ray mass recon¬ 
struction, including the effects of non-thermal pressure. The 
solid lines are the hydrostatic mass recovered from the X-ray 
measurements using Equation while the dashed lines are 
the hydrostatic mass plus a non-thermal component from 
Equation The shaded area is the 68% confidence limit 
allowed region from the weak lensing measurements, fit with 
an NEW profile. The vertical dotted line is the upper bound 
of our X-ray data; to the right of this line we have extrap¬ 
olated to illustrate the possible divergence of the mass esti¬ 
mates with and without significant non-thermal pressure. 

At all scales in Eigure the thermal pressure profile 
(solid line) is consistent with the shaded region, showing that 
the mass profiles estimated by the X-rays and lensing mass 
are consistent. This suggests that hydrostatic equilibrium is 
an acceptable approximation for our stacked profiles, given 
the error in our lensing measurements. 

We also see in Eigure that the thermal pressure profile 
with a non-thermal component (dashed line) enhances the 
hydrodynamical mass by 20% (10%) in the T < 2.5keV (T > 
2.5keV) cluster bin, but is still seen to be consistent with 
our lensing measurements. This shows that the non-thermal 
pressure expected from simulations falls within our present 
observed errors’, if present it acts in the opposite sense to 
chameleon gravity, reducing the detectable signal. 

With future X-ray measurements we will be able to fit 
out to a larger distance, allowing us to better constrain the 
effect of non-thermal pressure, which would be most promi¬ 
nent at large radii. We also note that our weak lensing pro¬ 
files have lower signal to noise than the X-ray profiles, how¬ 
ever with future lensing surveys we will be able to more 
accurately constrain these profiles also allowing us to bet¬ 
ter characterise not only chameleon gravity but non-thermal 
pressure too. 


6 CONCLUSIONS 

We have investigated the constraining power of stacked 
galaxy cluster profiles for testing chameleon gravity. We have 
examined 58 X-ray selected galaxy clusters, which have both 
good quality weak lensing data from CEHTlenS and X-ray 
data from XCS. After binning our clusters by X-ray tem¬ 
perature, we have generated weak lensing profiles and X- 
ray surface brightness profiles. Chameleon gravity predicts 
an additional pressure existing within clusters, which causes 
their gas component to become more compressed than GR 
gravity predicts. We have therefore investigated this phe- 
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Figure 5. Mass profile from the T < 2.5keV (T > 2.5keV) cluster bin in blue (red). The shaded area is the one-sigma allowed region 
from the weak lensing measurement and the solid line is the thermal mass reconstructed from the X-rays. The dashed line shows the 
thermal mass with an additional non-thermal component as discussed in the text. The vertical line is the upper extent of our X-ray data; 
to its right we have extrapolated the X-ray data. 


Scale 

Scale 

loglol/Rol 

Solar System 
(Hu & Sawicki|2007|) 

pc 

-6 

Dwarf Galaxies 
(Jain et al.|20I3|) 

kpc 

-6.3 

Coma cluster 
(Terukina et al.|20I4|) 

Mpc 

-4.2 

Cluster abundance 

Mpc 

—4.6 (n = I) 

(Cataneo et al.|20I4|) 

-3.5 (n = 3) 

Cluster stack 

Mpc 

-4.2 {n = I) 

(This Work) 

-3.7 {n = 3) 

CMB 

(Raveri et al.|20I4|) 

Gpc 

-3.0 


Table 1. Comparison of the constraints on log]^Q|/Ro|- 


nomena by comparing the X-ray profile with the weak lens¬ 
ing profile, which is unaffected by the fifth force. Using 
a multi-parameter MCMC analysis we have obtained con¬ 
straints on the common chameleon parameters [3 and 0oo, 
which in turn lead to constraints for |/ro|, a parameter char¬ 
actering f{R) theories. 


We find our results are competitive with other cosmo¬ 
logical constraints on chameleon models. In particular, our 
constraints are an order of magnitude stronger than those 



and |/ro|< 3.1xl0“^forn = 3 compared with our measure¬ 
ment of |/ro|< 2 X 10“^, all at the 95% CL. A comparison 
of these constraints is shown in Table [T] 

We examined the assumption of hydrostatic equilibrium 
by comparing the masses inferred from the X-ray observa¬ 
tions with weak lensing and found them to be consistent. 
Deviations from hydrostatic equilibrium would cause a dis¬ 
parity between the weak lensing and X-rays with the oppo¬ 


site sign to that from the chameleon effect. We modelled a 
non-thermal pressure X-ray component, and given current 
observational errors found this to be a subdominant effect 
on our constraints. 

As we are interested in the shape of the respective pro¬ 
files, the absolute mass of the stacked cluster, measured 
through both weak lensing and X-rays, is a nuisance pa¬ 
rameter which we have marginalised over. We therefore are 
not sensitive to the relative biases between these two tech¬ 
niques, such as reported in von der Linden et al. (2014) and 
(| 2 ^. 


Hoekstra et al. 


For the next generation of constraints via this method, 
we will need detailed modified-gravity hydrodynamic- 
simulations. These will allow us to check a range of as¬ 
sumptions used in this analysis such as hydrostaticity, non- 
thermal pressure, gas clumping in the cluster outskirts, 
spherical symmetry and the reliability of the NFW profile. 

We find our constraint on |/ro | to be consistent with the 
literature, and competitive at these cosmic scales and red- 
shifts. We have therefore demonstrated that it is possible to 
constrain chameleon gravity using stacked galaxy clusters; 
with the advent of wide area lensing surveys promising a 
much larger area, such as the Dark Energy Survey (DES, 
The Dark Energy Survey Collaboration|20Q5 ), the Kilo De¬ 


gree Survey (KIDS, de Jong et al. 2013), Euclid (Laureijs 
et al. 2011 1) and the Large Synoptic Survey Telescope (LSST, 


LSST Dark Energy Science Collaboration]|20I2 ), it will be¬ 


come possible to use stacks containing many more clusters 
to beat down systematics and obtain stronger constraints. 
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APPENDIX A: GOODNESS OF FIT 


To characterise the goodness of fit of our profiles we adopt 
the following statistic 

x" {tI , nl, b\, r}, Mloo , c’, Tq" , nf,', bj', r}', . (Al) 

II /O / \ 12, II 2 , 12, II 2 

C , P2, 000,2) — XWL + XWL + XSB + XSB 5 


where we adopt the notation I, H to indicate the temperature 
bins T < 2.5, T > 2.5 respectively, and 


I 2 

XwL — 2_^ 


(7(ri,) - 


(a’7,' 


,obs,I\ 2 


(A2) 


II 2 

XWL — 2 _^ 


(7(4',i) - 7U’")" 


(c^7f 


,obs,II\ 2 


(A3) 


xIb = X(^B(ri,0 - Stf)C-^{SB{r\,i) - S^Y), (A4) 

Xsb" = ^{SB{rlO-SY''')CrJiSBYL)-SYY- (A5) 

In the weak leasing case we approximate the covariance ma¬ 
trix as diagonal; we find strong leading diagonals for the 
measured correlation matrices. For the surface brightness 
fits we minimise over the full covariance matrix due to the 
covariances that exist between bins; here C is the error co- 
variance matrix. Then 7 (r^,i) is the value of the lensing 
model at a distance r± from the clusters’ centre; likewise 
SB(r±^i) is the value of the surface brightness model at a 
distance r± from the clusters centre. are the ob¬ 

served shear profile and surface brightness profile respec¬ 
tively, while is the observed error on the shear profile. 


APPENDIX B: SOURCE LIST 


XGS Name 

z 

Flag 

XMMXGS J020045.8-064229.2 

0.36 

0 

XMMXCS J020119.0-064954.6 

0.33 

0 

XMMXCS J020232.1-073343.8 

0.55 

1 

XMMXGS J020334.3-055049.5 


2 

XMMXCS J020359.1-055031.6 


3 

XMMXCS J020405.2-050142.5 


2 

XMMXCS J020428.5-070221.6 


2 

XMMXCS J020432.7-064449.4 


2 

XMMXCS J020514.7-045640.0 

0.29 

0 

XMMXCS J020611.4-061129.2 

0.88 

1 


Table Bl: Sample of the extended X-ray sources in 
GFHTFenS footprint. The XGS name and position are listed 
for all clusters. Redshifts are provided where available. The 
clusters forming the sample used throughout this work have 
a flag of 0 in the T < 2.5keV bin and a flag of 1 in the 
T > 2.5keV bin. A flag of 2 denotes the source was dis¬ 
counted for having no measured redshift. A flag of 3 denotes 
the source was discounted for having no measured X-ray 
temperature. The full version of this table is provided via 
the online edition of the article. An excerpt is provided above 
to illustrate form and content. 
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APPENDIX C: CLUSTER IMAGES 



Figure Cl: A selection of optically confirmed clusters as imaged by CFHTLenS. False colour composite images are 3' x 3b 
From left to right and top to bottom, the compilation shows the clusters: XMMXCS J020119.0-064954.6 at z=0.33; 
XMMXCS J021226.8-053734.6 at z=0.31; XMMXCS J021527.9-053319.2 at z=0.28; XMMXCS J021843.7-053257.7 at 
z=0.40; XMMXCS J022433.8-041433.7 at z=0.39; and XMMXCS J023142.2-045253.1 at z=0.21. These clusters are 
included in our sample, flagged either with a 0 or 1 in Table |B1| The remaining clusters in our compilation have no 
measured redshift or temperature and are flagged with a 2 or 3 in Table |B1[ Continuing onwards these clusters are: 
XMMXCSJ021517.1-0.60432.8; XMMXCSJ022359.2-083543.4; and XMMXCSJ141446.9+544709.1. 
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APPENDIX D: 2D CONTOURS 







E 


E 






I 


* 




m 


E 


E 


1 




« 


HBHH'J J 









m 





i 

i 


11 

■ 



a 

m 

01 


i 

i 

i 


■ 

■ 



^ ^ ^ cOr^cC^ cC^ C.y c> fN> 


„ II 

Uq 


c“ 


/?2 


<f>oc 


Figure Dl: The 95% (dark grey region) and the 99% CL (mid grey region) 2D marginalised contours for the 14 model 


parameters Tq [keV], Uq [10 


b\. 


[Mpc], M 200 [IO^'^Mq], c',To' [keV], n” [10 ^cm ®], b”,r“ [Mpc], M 2 I 0 


[lO^'^M©], c^,/32 , 000,2 used in our MCMC analysis. The rightmost plots show the ID likelihood distributions. 
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